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Abstract 

We present a study of the reaction e + e~ —>ir + ir~ tt° at the <p peak, W = M{(f)) 
= 1019.4 MeV, observed with the KLOE detector at DA<J>NE. The reaction is domi- 
nated by (f) production and decay, e + e~—Kp— >7r + 7r _ 7r . From a fit to the Dalitz plot 
density distribution we obtain the /3-meson parameters for its three charge states. 
We also find the relative amplitudes for (p—tpir and (f)—>Tr + ir~ ir and the cross section 
for e + e~ — ► lutt° with uj — > it + it~ . 

Key words: e + e~ collisions, <f> — > 7r + 7r~7r , p meson masses and widths 
PACS: 13.65.+i, 14.40.Cs 



The decay of the <fi meson to 7r + 7r~7r°, with a branching ratio (BR) of ~15.5%, 
is dominated by the pn intermediate states [1] p + n~, p~n + , and p°7r° with 
equal amplitudes. We use data from e + e - ^7r + 7r _ 7r° to determine the masses 
and widths of the three charge states of the p-meson. CPT invariance requires 
equality of the masses and widths of p + and p~ , while possible mass or width 
differences between p° and p^ 1 are related to isospin-violating electromagnetic 
effects. Additional contributions to e + e~— >7r + 7r _ 7r° are the so called "direct 
term", 0^7r + 7r~7T°[2], and e + e~ — > um°, uj — > ty + ti^. 
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We define the variables x = T + — T~ and y = T°, where T + ' ~' are the 
kinetic energies of the three pions in the center of mass system (CM). The 
Dalitz plot density distribution D(x,y), in terms of the three amplitudes is 
given by: 

D(x,y) oc \pl x pl\ 2 \A p7T + A dir + A un \ 2 (1) 

where p± are the 7r ± momenta in the CM and A p7T , A^n and A W7r are the 
amplitudes described above. 

Some 2 million ir + 7i^7i° events were collected with the KLOE detector in e + e~ 
collisions at the Frascati 0-factory DA$NE, in the fall of 2000. DA$NE was 
run at a CM energy W = M(0) and the integrated luminosity was L=16 pb^ 1 . 
In the following we present an analysis of the Dalitz plot density distribution 
and new determinations of T(p + ~°), M(p +_0 ), A p7T , Aa T , A^ and cr(e + e~ — > 
uti ,uj -> 7T+7T-) at W = 1019.4 MeV. 

The KLOE detector consists of a large volume drift chamber [3] (3.3 m length 
and 2 m radius), operated with a 90% helium-10% isobutane gas mixture, 
and a sampling electromagnetic calorimeter [4] made of lead and scintillating 
fibres. The whole detector is surrounded by a superconducting coil produc- 
ing a solenoidal field 5=0.52 T. The drift chamber momentum resolution is 
cr(p_i_)/p_i_~0.4%. The calorimeter determines photon impact points to an ac- 
curacy of 1 cm / y / £'(GeV) in the direction along the fibres and of 1 cm in 
the transverse direction. Photon energies and arrival times are measured with 
resolutions of a{E)/E = 5.7%/y / £(GeV) and a(t) = 57ps/y / E(GeV) ©50 ps. 
The solid angle coverage of the calorimeter is 98% of 47r; however particles 
from the interaction point with a polar angle with respect to the beam axis 
9 < 15° are shadowed by the low-/3 insertion quadrupoles. The trigger [5] is 
based on the coincidence of at least two local energy deposits in the calorime- 
ter, above a threshold that ranges between 50 and 150 MeV. In order to reduce 
the trigger rate due to cosmic rays crossing the detector, events with a large 
energy release in the outermost calorimeter planes are vetoed. 

The DA$NE beams collide with a crossing angle of 25 mrad, so the is pro- 
duced with a momentum in the horizontal plane [p^| ~ 12.5 MeV. 7r + 7r~7r° 
events are selected by asking for two non-collinear tracks with opposite sign 
of curvature and polar angle 9 > 40° which intersect the interaction region. 
The acollinearity cut (A# < 175°) removes e + e~7 events without incurring an 
acceptance loss for the signal. We then compute the missing mass, 
M miss = J (E^ — En+ — E v -) 2 — — p v + — p n - 1 2 where E and p are labora- 
tory energies and momenta. M m i SS is required to be within 20 MeV of the 7r° 
mass. This requirement corresponds to an effective energy cut of < 20 MeV on 
the total energy radiated because of initial state radiation (ISR). Two photons 
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M„ (MeV) cosi5„ 

Fig. 1. Distributions of M m j ss (top), M 77 , and cos# 77 (bottom left and right) for a 
sample of selected events. The rms widths of the M m j ss and M 77 distributions are 
5.5 MeV and 17 MeV, respectively. The solid lines are gaussian fits. 

in the calorimeter are also required. A photon is defined as an energy deposit 
larger than 10 MeV with 21° < 9 < 159° and an arrival time compatible with a 
particle travelling at the speed of light, within 5a(t). The two-photon opening 
angle in the ir° rest frame must satisfy cos^ < —0.98. 

Fig. 1 shows the distributions of the missing mass M miss , of the 77 invariant 
mass, and of cos# 77 for a sample of selected events. Due to the large cross- 
section 2 for this final state with respect to other processes (cx^, x BR(0 — > 
7r + 7r~7T°) = 460 nb) and to the clean signature, the background to this process 
after the selection described is < 10 -5 . The Dalitz plot variables x and y are 
evaluated using the measured momenta of the charged pions, boosted to the 
center of mass system: x = E*^ — E*_ and y = E^ — E+ — E*_ — M n o = T n o . E^ 
and Prf, are measured run by run using Bhabha scattering events. ISR lowers 
the mean 7r + 7r~7r° total energy by ~130 keV. This value is used in the analysis 
with negligible effect on the results. The resolution on x and y is about 1 MeV 
over the full kinematical range. 

The Dalitz plot density distribution is shown in Fig. 2. In the plot the num- 
ber of events corrected for the efficiency is shown divided by \p+ x pl\ 2 . 
Three bands corresponding to the three p states are clearly evident. The two- 



2 Here and in the following we consider visible cross-sections, not corrected for the 
effect of the radiative corrections. 
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Fig. 2. Distribution of the number of events corrected for the efficiency and divided 
by \pl x p*\ 2 . The grey scale is in arbitrary units. The plot contains 1.98 millions 
events in 1874 bins 8.75x8.75 MeV 2 each. Three broad bands corresponding to the 
three p states are indicated. The kinematical boundary is also shown. 

dimensional distribution is plotted in 8.75x8.75 MeV 2 bins. There are 1874 
bins within the kinematic boundary. The bin width is larger than the x and y 
resolution, but is small compared to the density variations of the Dalitz plot 
as can be seen in the x and y projections shown in Fig. 3. Smearing effects 
due to the resolution are negligible. 



Trigger and selection efficiencies have been evaluated as functions of x and y. 
A full Monte Carlo simulation of the detector has been used with corrections 
based on control samples of data. Corrections to the detection efficiency for 
low energy photons have been obtained using e + e~7 events while analysis of 
tracking efficiency shows that the ratio of data to Monte Carlo efficiency is very 
close to 1 and uniform in the relevant momentum range. The trigger cosmic 
ray veto rejects ~5% of the 7r + 7r~7r° events. The amount and the distribution 
of the rejected events in the Dalitz plot has been evaluated using downscaled 
samples of vetoed events. Fig. 4 shows the efficiency bin by bin obtained from 
a sample of 5 x 10 6 Monte Carlo events. It can be noticed that the overall 
efficiency ranges between 20% and 30% within the Dalitz plot. The efficiency 
is dominated by acceptance cuts on charged tracks and photons. In particular, 
the low efficiency in the upper corners is due to the fact that low momentum 
pions do not reach the drift chamber. The uncertainty on the efficiency is 
dominated by Monte Carlo statistics and is on the order of 1% in each bin. 

The calibration of the momentum scale is checked using the measured missing 
mass M miss . The central value of M miss differs from M n o by less than 0.1 MeV. 
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Fig. 3. Distributions in x and y of the number of events. The u peak is not visible 
because of the smallness of the signal. 

The rms fluctuation of the missing mass across the entire Dalitz plot is 0.25 
MeV. The latter value is used in the estimates of systematic uncertainties for 
masses and widths. 



The Dalitz plot distribution has been fitted according to Eq. 1. The pir term 
includes the line-shape of the three charge states of the p [6,7,8] : 



.4 



p7T 



Ml 



ql-Ml + iq k T {k) {ql) 



(2) 



where k = +, — , is the p-meson charge and F^iOk) * s defined as: 



r. 



' ftr(gfc) 

M M l) 



'Ml 

, 4 , 



(3) 



In the above, q\ is the invariant mass of the appropriate pion pair, 7rV J with 
i, j, k = +, — , 0; Mk is the p mass; the width; and p n is the pion momentum 
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Fig. 4. Overall trigger and selection efficiency including the applied corrections as 
a function of the position in the Dalitz plot. 

in the p center of mass. The direct amplitude does not in general have the 
same phase as A pn and is taken as A^ iT = a^e*^ . This term can be interpreted 
as the sum of p-like resonances of higher mass like <p p'(1415)7r, which is not 
distinguishable from a constant term multiplied by x xp*_\ 2 . The uir term 
is restricted to a horizontal band around y = 83.7 MeV, the kinetic energy of 
the 7T° in the ujtt final state. The amplitude is given by: 

M 2 

A - n pi4>w _Zj£ (A\ 

Qo~M 2 + tqoT^ 



with the uj parameters fixed to PDG values [9]. a p , and are three dimen- 
sionless parameters. 

We perform a \ 2 n t °f the described density function to the data binned as 
described. The errors used are from data statistics and Monte Carlo statistics 
for efficiency calculations. Three fits have been performed: a) a fit assuming 
CPT and isospin invariance, i.e. M p o = M p + = M p -, T p o = T p + = T p -; b) a 
fit assuming only CPT invariance i.e. M p + = M p -, T p + = T p -; and finally c) 
without limitations on masses and widths. Moduli and phases of direct and urn 
terms and a p are also free parameters. There are 7 free parameters for fit a), 9 
for fit b) and 11 for fit c). The number of degrees of freedom is 1874—7=1867 
for fit a), 1865 for fit b) and 1863 for fit c). The results of the fits are shown 
in Tab. 1 and a one-dimensional comparison of the data with the fit function 
(fit (c)) is shown in Fig. 5. 

The uncertainties on the parameters are given as the sum of fit uncertainties 
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Fig. 5. Comparison between efficiency corrected data (points, N/e is number of 
events per bin divided by the bin efficiency) and fitted function (histogram) as a 
function of the bin number. The structure observed is due to the y distributions for 
x slices. 

(due to event statistics and efficiencies) and systematic uncertainties coming 
from two sources: (a) the stability of the fit with respect to changes of the selec- 
tion cuts, and (b) the absolute momentum calibration. Source (a) is dominant, 
and is evaluated by repeating the fit for data samples obtained using different 
values of the cuts and taking the rms variation of the parameters. Effects due 
to the modelling of ISR are automatically taken into account by the variation 
of the cut on the missing mass. 

Fit (a) converges to an acceptable \ 2 probability, showing that the experi- 
mental distribution is consistent with CPT and isospin invariance. Mass and 
width differences between states of charge obtained by fit (b) and fit (c) are 
summarized in Tab. 2 where M p ± 



(M p+ +M p -)/2 and T p ± = (r p+ + r p -)/2. 
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Table 1 

Fit results. Masses and widths are in MeV. The amplitudes are scaled in such a 
way that a p =l. 4>d and <f> w are in radians. % 2 values and probabilities of the fits are 
given in the first row. 



parameter 


fit (a) 


fit(b) 


fit(c) 


A LrVA )\ 


1939 [12%1 


1914 [21%1 


1902 [26%1 






775.9 ± 0.5 ± 0.5 


775.9 ± 0.6 ± 0.5 


M + 


775.8 ± 0.5 ± 0.3 


775.5 ± 0.5 ± 0.4 


776.3 ± 0.6 ± 0.7 


M n - 






774.8 ± 0.6 ± 0.4 


T„o 




147.3 ± 1.5 ± 0.7 


147.4 ± 1.5 ± 0.7 




143.9 ±1.3 ±1.1 


143.7 ±1.3 ±1.2 


144.7 ± 1.4 ± 1.2 








142.9 ± 1.3 ± 1.4 


ad 


0.78 ±0.09 ±0.13 


0.72 ± 0.09 ± 0.05 


0.71 ±0.09 ±0.05 


(t>d 


-2.47 ±0.08 ±0.08 


-2.43 ±0.10 ±0.08 


-2.43 ±0.10 ±0.08 


x 10 3 


7.1 ±0.6 ±0.8 


9.0 ±0.7 ±0.7 


9.0 ±0.7 ±0.3 


<t>LU 


-0.22 ±0.11 ±0.04 


-0.10 ±0.10 ±0.05 


-0.10 ±0.10 ±0.07 



The p masses are significantly larger and the widths smaller than the PDG 
averages (M p = 771.1 ± 0.9 MeV and T p = 149.2 ± 0.7 MeV, see [9]) but are 
close to recent results [10,11,12] for the p° mass and width. 

Table 2 

Mass and width differences (in MeV) between p- mesons, from fit (b) and fit (c). 
Errors include the correlations between parameters. 





-M p ± 


0.4 ±0.7 ±0.6 


M p+ 


-M p - 


1.5 ±0.8 ±0.7 






3.6 ±1.8 ±1.7 






1.8 ±2.0 ±0.5 



We observe significant contributions from the direct term and from the uir 

term: the coefficients and are significantly different from 0. If we define 

the weight of each contribution a (a = pir, dir and ujtt) as 

I a — I dxdy\A a \ 2 / J dxdy\A tot \ 2 , where the amplitudes A a are taken from the 

results of fit c) and A tot is the sum of three amplitudes, we get the following 

weights: 

l lm = 0.937 (5) 
J dir = 8.5 x 10~ 3 (6) 
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1^ = 2.0 x 10~ 4 (7) 

The sum of the three weights is not equal to one, due to a sizeable interference 
term between A pv and A^ ir accounting for about 6% of the total. 

From I wn , if we neglect interference, we obtain the visible cross section for the 
non resonant process e + e~ — * ujtt with u — > 7r + 7r~ at W — 1019.4 MeV. 

a(e + e- -> un -> tt+tTtt ) = I U7r x a(e + e~ -> Tr+vr^Tr ) = 92 ± 15 pb(8) 

The ratio of this cross-section to the value <j(e + e _ — > ujtt — > 7r°7r°7) = 0.46 ± 
0.01 ± 0.03 nb obtained by KLOE in the 5 photons analysis [13] gives: 

BR(u -> TT+TT") 

which compares well with the value 1.7/8. 7=0. 20±0. 03 from [9]. 
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